Reduction of Two-Dimensional Dilute Ising Spin Glasses 
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The recently proposed reduction method is applied to the Edwards-Anderson model on bond- 
diluted square lattices. This allows, in combination with a graph-theoretical matching algorithm, 
to calculate numerically exact ground states of large systems. Low-temperature domain-wall ex- 
citations are studied to determine the stiffness exponent yi- A value of yi = — 0.28f(3) is found, 
consistent with previous results obtained on undiluted lattices. This comparison demonstrates the 
validity of the reduction method for bond-diluted spin systems and provides strong support for 
similar studies proclaiming accurate results for stiffness exponents in dimensions d = 3, . . . , 7. 
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I. INTRODUCTION 

Despite more than two decades of intensive research, 
many properties of spin glasses 0, 0, 01 3 , especially in 
finite dimensions, are still not well understood. The most 
simple model is the Edwards- Anderson model (EA) Q , 



H 



(i) 



<i,3> 



with Ising spins Xi = ±1 arranged on a finite-dimensional 
lattice with nearest-neighbor bonds Ji_j, randomly drawn 
from a distribution P{J) of zero mean and unit variance. 

For two-dimensional Ising spin glasses it is now widely 
accepted that no ordered phase for finite temperatures 
exists UlSHSli, wnne spin glasses order at low tem- 
peratures in higher dimensions |lfj . ITTl Ir3 , Il3l | . This can 
be seen e.g. by studying the stiffness exponent y (often 
labeled 6) which is in many respects one of the most fun- 
damental quantities to characterize the low-temperature 
state of disordered systems. This exponent provides an 
insight into the effect of low-energy excitations of such a 
system 0,0. A recent study suggested the importance 
of this exponent for the scaling corrections of many ob- 
servables in the low-temperature regime [T^ | , and it is an 
essential ingredient to understand the true nature of the 
energy landscape of finite-dimensional glasses [l7tllcill9| . 

To illustrate the meaning of the stiffness exponent, one 
my consider an ordinary Ising ferromagnet of size L d hav- 
ing bonds J = +1, which is well-ordered at T — for 
d > 1, with periodic boundary conditions. If we make 
the boundary along one spatial direction anti-periodic, 
the system would form an interface of violated bonds 
between mis-aligned spins, which would raise the en- 
ergy of the system by AE ~ This "defect" - 
energy AE provides a measure for the energetic cost 
of growing a domain of overturned spins, which in a 



ferromagnet simply scales with the surface of the do- 
main. In a disordered system, say, a spin glass with 
an equal mix of J — ±1 couplings, the interface of 
such a growing domain can take advantage of already- 
frustrated bonds to grow at a reduced or even decreasing 
cost. To wit, we measure the probability distribution 
P(AE) of the interface energies induced by perturba- 
tions at the boundary of size L, for which the typical 
range a{AE) = V< AE 2 > — < AE > 2 of the defect- 
energy may scale like 



a(AE) ~ L v . 



(2) 
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From the above consideration it is clear that y < d — 1 , 
and a bound of y < (d — l)/2 has been proposed for 
the Edwards- Anderson model generally Particular 
ground states of systems with y < would be unstable 
or only marginally stable with respect to spontaneous 
fluctuations, whether induced thermally or structurally. 
These fluctuations could grow at no cost, like in the case 
of the one-dimensional ferromagnet where y = d— 1 = 0. 
Such a system does not manage to attain an ordered state 
for any finite temperature. Conversely, a positive sign 
for y at T — indicates a finite-temperature transition 
into an ordered regime while its value is a measure of the 
stability of the ordered state. It is generally believed that 
the EA possesses a glassy low-temperature regime, i. e. 

t> o, for d> 3\&MMMmMMM\mWl 
, while such a phase is absent, and y < 0, for d < 2 
123 . l2q . A value of y — marks the lower critical 
dimension. 

The stiffness exponent provides a measure of the ef- 
fect of excitations on a spin system around ground state 
configurations, induced by low-temperature fluctuations. 
It is computationally convenient to induce such excita- 
tions by perturbing the system of size n = L d at one 
boundary and measuring the response for increasing sys- 
tem size L. The square lattices considered here have 
one open and one periodic boundary, and we determine 
the energy difference AE between the ground states of 
a given bond-configuration, once with periodic and once 
with anti-periodic boundary conditions (P- AP method) . 
Anti-periodic boundary conditions are obtained by re- 
versing the sign of a strip of bonds along the periodic 
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boundary. Using a symmetric bond distribution P(J), 
AE will be also symmetrically distributed, but with a 
variance a of these excitation that may scale with L ac- 
cording to Eq. (J2J. 

Until recently, the consensus of the results for d = 3 
ranged from y 3 m 0.19 to « 0.27 [H GJ, E III SJ El 
l23l |3 H3 ) while there was only two results in d = 4, 
2/4 = 0.64(5) [11 and y 4 = 0.82(6) H3- In Refs. [MIH, 
it was proposed to study the EA in Eq. JIJ on bond- 
diluted lattices to obtain more accurate scaling behavior 
for low-temperature excitations. One can remove itera- 
tively low-connected spins from the lattice and alter the 
interactions, i.e. reduce the system, in such a way that the 
ground-state energy of the reduced system is the same as 
the original system. In this way often much larger lattice 
sizes L can be simulated compared to undiluted ones and, 
in combination with finite-size scaling, enhanced scaling 
regimes are achieved. In this manner, improved or en- 
tirely new values for the T = stiffness exponents in di- 
mensions d = 3 to 7 were computed for lattices with dis- 
crete bonds, ±J, resulting in 2/3 = 0.24(1), 2/4 = 0.61(2), 
7/5 = 0.88(5), y 6 = 1.1(1), and y 7 = 1.24(5). 

The novelty of the procedure used in Refs. |2^, |3(j 
makes it difficult to assess the validity and the accuracy 
of the approach, since few data of comparable accuracy 
exists for the results presented there. The scaling Ansatz 
used is based on various reasonable but untested assump- 
tions, for instance that the stiffness exponent does not de- 
pend on the dilution. Therefore, this approach has been 
discussed in the context of the Migdal-Kadanoff approxi- 
mation [inj to justify the scaling Ansatz. Here, we apply 
this procedure to the EA in d = 2, which has been studied 
extensively in recent years. These studies have found the 
value of the stiffness exponent y to be 2/2 = —0.281(2) 
y 2 = 0.287(4) |H, or y 2 = -0.282(2) [l3|. We 
find that the result obtained here on diluted lattices, 
1J2 = —0.281(3), compares well with those earlier re- 
sults. This does not add any accuracy to the value of 
2/2, and we find anew |3l| that diluted lattices with con- 
tinuous (Gaussian) bonds are beset with more complex 
scaling behavior |32( as well as more extensive scaling 
corrections. Much more important, though, the present 
study indicates the correctness of the reduction approach, 
hence the validity of the results obtained for larger di- 
mensions. In particular, the stiffness exponent does not 
depend on the dilution, i.e. it is universal, even when 
2/<0. 

In the next section, we describe the algorithm we ap- 
plied to reduce and evaluate large instances of dilute, pla- 
nar lattices and outline the ground-state algorithm. In 
Sec. IIVI we discuss the results of our numerical studies, 
followed by some conclusions in Sec. 



II. ALGORITHMS 

In this section we explain the algorithms used to calcu- 
late exact ground states of diluted two-dimensional Ising 



spin glasses. Our approach consists of two steps, each 
previously introduced in Refs. |3(j and 0, which we 
will review in the following. First, the systems are re- 
duced, i.e. low-connected spins are iteratively removed, 
while altering the remaining interactions such that the 
ground-state energy is not affected. After the reduction 
is finished, the ground state of the remaining system is 
calculated exactly using a matching algorithm. 

Absent a true glassy state in d = 2, it is not too sur- 
prising that computationally efficient ground-state algo- 
rithms exist [33I l34l l35| which exhibit a running time 
growing only polynomially with system size. This allows 
to measure y with great accuracy. For d > 3, where a 
true glassy state exists at low temperatures, no computa- 
tionally efficient methods are known to determine ground 
states exactly. The ground-state calculation belongs j3|J 
to the class of NP-hard problems [37|, where all exist- 
ing algorithms exhibit an exponentially growing running 
time with size. Instead, heuristic optimization methods 
pM l39j | are used which typically are believed to approxi- 
mate ground states for lattices with up to n « 10 3 spins 
(or L < 14 in d = 3) with some confidence. Any inac- 
curacy in the determination of such ground states gets 
further aggravated by way of the subtraction leading to 
AE, suggesting that the scaling regime in d = 3 extends 
at most up to L sa 10, and even less in d > 3. 

In light of those difficulties, it might come as a sur- 
prise that the study of the EA on diluted lattices would 
possibly improve matters. After all, dilution eases the 
constraintness of the spin configuration, leading to less 
frustration, and locally to a less glassy state at low tem- 
peratures. Consequently, the length scale beyond which 
frustration effects local spin arrangements should be ex- 
tended for increasing dilution, leading to persistent scal- 
ing corrections before an asymptotic scaling regime can 
be obtained at much larger system sizes. Thus, any gain 
in obtainable system size provided by the dilution should 
only marginally effect any useful scaling regime. Yet, 
the numerical results using a ± J bond-distribution prove 
otherwiseRefs. [2!3,|3(j: While scaling corrections worsen 
as expected at too small bond-densities, they are sig- 
nificantly suppressed at intermediate densities even com- 
pared to the undiluted case. The origin of those reduced 
scaling corrections at intermediate densities has been in- 
vestigated in Ref. [3]]]. Additionally, collapsing all data 
from various bond fractions p with a scaling Ansatz ex- 
tends scaling even further. 



A. Reduction Method 

To exploit the advantages of spin glasses on a bond- 
diluted lattice, we can often reduce the number of rele- 
vant degrees of freedom substantially before a call to an 
optimization algorithm becomes necessary. Such a re- 
duction, in particular of low-connected spins, leads to a 
smaller, compact remainder graph, bare of trivially fluc- 
tuating variables, which is easier to optimize. Here, we 
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focus exclusively on the reduction rules for the energy at 
T = 0; a subset of these also permit the exact determina- 
tion of the entropy and overlap ■ These rules apply to 
general Ising spin glass Hamiltonians as in Eq. with 
any bond distribution P(J), discrete or continuous, on 
arbitrary sparse graphs. 

The reductions effect both spins and bonds, eliminat- 
ing recursively all zero-, one-, two- ; and three-connected 
spins. From previous applications [2!|,[3(|, we have sup- 
plemented these rules with one that is not topological but 
concerns bond values directly, which is especially effec- 
tive for Gaussian bond distributions. These operations 
eliminate and add terms to the expression for the Hamil- 
tonian in Eq. ft), but leave it form-invariant. Offsets in 
the energy along the way are accounted for by a variable 
H 0l which is exact for a ground-state configuration. 

Rule I: An isolated spin can be ignored entirely. 

Rule II: A one-connected spin i can be eliminated, 
since its state can always be chosen in accordance with its 
neighboring spin j to satisfy the bond JJj ,-. For its ener- 
getically most favorable state we adjust H := H — \Ji.j | 
and eliminate the term — Jy Xi Xj from H. 

Rule III: A double bond, J$ and J- 2 . , between two 
vertices i and j can be combined to a single bond by 
setting Jij — J>j + jf9 or be eliminated entirely, if the 
resulting bond vanishes. This operation is very useful to 
lower the connectivity of i and j at least by one. 

Rule IV: For a two-connected spin i, rewrite in Eq. Q 

Xi(Ji,lXx + Ji, 2 X 2 ) < | Ji,lX\ + Ji^2X2\ 

= Ji, 2 x x x 2 +AH, (3) 

where 

Jl,2 = -^{\Ji,l + J'i,2\ — \Ji,l — Ji,2\) , 

AH = l(\Ji,l + Ji, 2 \ + \Ji,l-Ji,2\), (4) 

leaving the graph with a new bond Ji f2 between spin 1 
and 2, and acquiring an offset H Q := H Q — AH. 

Rule V: A three-connected spin i can be reduced via a 
"star-triangle" relation, as depicted in Fig. ^ 

Xi X2 + Ji,3 

Xi x 3 

< I Ji,\X\ + Ji, 2 X 2 + Ji,3X 3 \ (5) 
= Ji >2 Xi X 2 + Ji, 3 Xi X 3 + J 2 ,3 X 2 X 3 + AH, 

where 

Ji, 2 = -A - B + C + D, Ji, 3 = A - B + C - D, 
J 2 , 3 = -A + B + C-D, AH = A + B + C + D, 

A = 4 \ Ji,l — Ji,2 + Ji,3\ i B = j \ JiA — Ji.2 — Ji,3\ , 
C = 4 \ J%,1 + Ji,2 + Ji,3\ i D = j | Jj ( l + Ji_2 — Ji, 3 \ ■ 

Rule VI: A spin i (of any connectivity) for which the 
absolute weight \Ji,j>\ of one bond to a spin j' is larger 




FIG. 1: "Star-triangle" relation to reduce a three-connected 
spin xo- The new bonds on the right are obtained in Eq. JSJ. 




FIG. 2: Illustration of Rule VI for "strong" bonds. Left, 
the local topology of a graph is shown for two spins, xo and 
xi, connected by a bond Jo,i (thick line). If Jo,i > (resp. 
Jo,i < 0) satisfies Eq. @, xo and xi must align (resp. anti- 
align) in the ground state and xo can be removed. Right, the 
remainder graph is show after the removal. The other bonds 
emanating from xo (dashed lines) are now directly connected 
to x\ (potentially with a sign change, if Jo,i < 0). This 
procedure may lead to a double bond (Rule III), if xi was 
already connected to a neighbor of xq before. 

than the absolute sum of all its other bond-weights to 
neighboring spins j ^ j', i. e. 

i^i>Ei j ^i> ( 6 ) 

bond Ji ji must be satisfied in any ground state. Then, 
spin i is determined in the ground state by spin j' and it 
as well as the bond can be eliminated accordingly, 
as depicted in Fig- El Here, we obtain H := Hq — \Jij>\- 
All other bonds connected to i are simply reconnected 
with j', but with reversed sign, if < 0. 

This procedure is costly, and hence best applied after 
the other rules are exhausted. But it can be highly ef- 
fective for very widely distributed bonds. In particular, 
since neighboring spins may reduce in connectivity and 
become susceptible to the previous rules again, further 
reductions may ensue, see Fig. 

The bounds in Eqs. IMKil) become exact when the re- 
maining graph takes on its ground state. Reducing even 
higher-connected spins would lead to new (hyper-)bonds 
between more than two spins, unlike Eq. ft). While such 
a reduction is possible and would eventually result in the 
complete evaluation of any lattice ground state, it would 
lead along the way to an exponential proliferation in the 
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number of (hyper-)bonds in the system. This fact is a 
reflection of the combinatorial complexity of the glassy 
state, which will be explored in Ref. |4fJ. 

After a recursive application of these rules, the origi- 
nal lattice graph is either completely reduced (which is 
almost always the case below or near p c ), in which case 
H Q provides the exact ground state energy already, or we 
are left with a highly reduced, compact graph in which 
no spin has less than four connections. We obtain the 
ground state of the reduced graph with an exact match- 
ing algorithm as used in Ref. [jj , which together with H a 
provides the ground-state energy of the original diluted 
lattice instance. 
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FIG. 3: 2d spin glass with all spins up (left, up spins not 
shown). Straight lines are ferromagnetic, jagged lines are 
anti-ferromagnetic bonds. The dotted lines connect frustrated 
plaquettes (crosses). The bonds crossed by the dotted lines 
are unsatisfied. In the right part the ground state with three 
spins pointing down (all other up) is shown, corresponding to 
a minimum number of unsatisfied bonds. 



B. Matching 

Let us now explain just the basic idea of the match- 
ing algorithm, for the details, see Refs. [33|,02 EBl- The 
method works for spin glasses which are planar graphs; 
this is the reason, why we apply periodic boundary con- 
ditions only in one direction. In the left part of Fig. [3| 
a small 2d system with open boundary conditions is 
shown. All spins are assumed to be "up" , hence all anti- 
ferromagnetic bonds are not satisfied. If one draws a dot- 
ted line perpendicular to all unsatisfied bonds, one ends 
up with the situation shown in the figure: all dotted lines 
start or end at frustrated plaquettes and each frustrated 
plaquette is connected to exactly one other frustrated 
plaquette. Each pair of plaquettes is then said to be 
matched. Now, one can consider the frustrated plaque- 
ttes as the vertices and all possible pairs of connections 
as the edges of a (dual) graph. The dotted lines are se- 
lected from the edges connecting the vertices and called a 
perfect matching, since all plaquettes are matched. One 
can assign weights to the edges in the dual graph, the 
weights are equal to the sum of the absolute values of 
the bonds crossed by the dotted lines. The weight A 
of the matching is defined as the sum of the weights of 
the edges contained in the matching. As we have seen, 



A measures the broken bonds, hence, the energy of the 
configuration is given by E — — Yluj) + 2 A. Note 
that this holds for any configuration of the spins, since a 
corresponding matching always exists. Although Fig. [31 
only shows a square lattice, a matching is always possible 
for any planar graph, such as the reduced, dilute lattices 
discussed here. 

Obtaining a ground state means minimizing the total 
weight of the broken bonds (see right panel of Fig. |3J), 
so one is looking for a minimum- weight perfect matching. 
This problem is solvable in polynomial time. The algo- 
rithms for minimum- weight perfect matchings [4lL |42| are 
among the most complicated algorithms for polynomial 
problems. Fortunately the LEDA library offers a very ef- 
ficient implementation |43| . which we have applied here. 



III. SCALING ANSATZ 

Clearly, there exists a lowest bond fraction p* , below 
which a glassy state is not possible. In particular, the 
lattice must exceed the bond-percolation threshold p c to 
exhibit any long-range correlated behavior and p* > p c 
must hold. It is expected that p* = p c for any continuous 
distribution, while for discrete distributions p* may be 
minutely larger thanp c [2^,|2^,|32. Accordingly, for two- 
dimensional bond-diluted lattices and a Gaussian bond 
distribution used here, we expect p* = p c = 1/2 [44| . 
Similarly, it is expected that the correlation length near 
the critical point scales as 



Up)~{p-p*)- v 



(7) 



with v* = v = 4/3, well-known from 2-dimensional per- 
colation. 

The introduction of the bond density p as new pa- 
rameter permits a finite-size-scaling Ansatz in the limit 
P — > p*. Combining the data for all L and p leads to a 
new variable x = L/£(p), which has the chance of ex- 
hibiting scaling over a wider regime than for L alone. As 
we have argued in Ref. [3l|], this Ansatz should take the 
form 



a(AE) LtP ~ £(p)» 



/(*), 



(8) 



as suggested by Refs. .32-, 45). In principle, the Ansatz 
requires L 1 and £(p) ^> I. The scaling function / 
was chosen as to approach a constant for L 3> £,(p)- 

Note that one basic assumption used here is that the 
exponent y does not depend on the bond density p. Since 
a percolating cluster sufficiently above the percolation 
transition is effectively compact for large L, one may ar- 
gue that asymptotic scaling properties of the spin glass 
should be uneffected by p. Reproducing the scaling of the 
undiluted lattice with the Ansatz in Eq. (jHJ with some 
accuracy would add support to this argument. 

To obtain the exponent yp in Eq. © directly, one con- 
siders the limit L ~ £(p) — > oo, i. e. x ~ I. Then, one 
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can show that at p = p c |31j, |4 



a(AE) 



L y r 



(9) 



We observe that due to the fractal nature of the per- 
colating cluster at p* = p c , no long-range order can 
be sustained and defects possess a vanishing interface. 
Thus, yp < and the defect energy vanishes. In con- 
trast, yp = for discrete bond distributions such as ± J, 
because p* > p c and the interface energy a(AE) be- 
comes L-independent. As it turns out, with yp = the 
scaling collapse greatly simplifies for ±J bonds, leaving 
continuous bonds with one extra exponent to account for. 
Hence, using Gaussian bonds, the accuracy obtained for 
the desired stiffness exponent y diminishes, as the study 
in Ref. |3J shows. 

In case of a finite-temperature glass transition with di- 
vergent energy scales (y > 0), universality provides us 
with the choice of the more convenient distribution, ±J, 
to compute the stiffness exponent [2|| |3(| • Apparently, 
this universality brakes down below the lower critical di- 
mension, d < di ss 2.5 |sl l46j . and a nontrivial value for y 
is only obtained for continuous bonds y|. Thus, we have 
to take the exponent yp in Eq. JSJ) into account. 

A scaling collapse is further complicated by the fact 
that the asymptotic regime of interest for the determina- 
tion of y, namely x >■ 1 or L 3> £(j>), is hard to access 
for p — > p* . Most data that reaches asymptotic scaling, 
i. e. x ^> 1, is typically obtained instead at intermediate 
values of p, sufficiently above p* to reach system sizes 
with L ^ but small enough to exploit the reduction 
rules from Sec. Ill Al As the analysis in Ref. [3l| sug- 
gest, in that regime the correlation length Eq. Q) may 
be too small to justify Eq. JSJ. Furthermore, it is clear- 
that Eq. Q is valid only close to p Cl leading possibly to 
wrong estimations for critical exponents obtained from 
finite-size scaling 0] • Also the finite-size corrections for 
the correlation length itself play probably an important 
role, which can bee seen from previous studies |4S| of 
two-dimensional percolation, where a significant change 
of the effective exponent v was observed when changing 
the system size from L = 2 to L = 10000. Similarly, un- 
known scaling corrections missing in the form of Eq. © 
are likely to arise. Yet, experience shows that a focus on 
data with L ^> £(p) for any £(p) at least provides a sat- 
isfactory collapse in the i» 1 regime with an accurate 
prediction for the stiffness exponent y. There, unlike y, 
the exponents v* and yp, and the scaling function f(x), 
which are more closely associate with the scaling window 
near p* , will not be accurately represented. In any fit of 
the data, their values are likely distorted to absorb the 
effect of unknown scaling corrections. 

For our data analysis, we will therefore apply cuts to 
eliminate data outside of x 3> 1, and fit the remaining 
data to the form 



<7(A£)l, p 



/(oo), 



(10) 



D 




0.01 r 



fixing = (p-p* 



with /(oo), p* , v* , yp, and y as 



FIG. 4: Logarithmic plot of gl at p = p* — 1/2. An asymp- 
totic power-law regime corresponding to Eq. @, different 
from those for p > p* , is reached quickly, and an asymptotic 
fit extrapolates to yp = —0.98(2). 



fitting parameters. Note that in this limit, v* and yp are 
not independent. Hence, we choose to fix v* — v = 4/3. 



IV. NUMERICAL EXPERIMENTS 

In our numerical simulation, we have generated a large 
number of instances of symmetric Gaussian bond disor- 
der on square lattices with open boundaries vertically, 
and periodic boundary conditions horizontally, as de- 
scribed in Ref. Q. But these instances are bond-diluted 
with a bond fraction of p > p* = p c — 1/2. On this bond- 
diluted spin glass, the reduction algorithm from Sec. Ill Al 
is applied to recursively remove as many spin variables 
as possible while exactly accounting for their contribu- 
tion to the ground-state energy. Since the original lattice 
was a planar graph, the reduction method preserves this 
property for the remainder graph. Hence, the matching 
algorithm discussed in Sec. Ill Bl can be applied to the 
remainder graphs here to determine their exact ground 
states in polynomial time. In this manner, we study the 
defect energy a(AE) both as a function of size L and 
bond density p. We consider systems of sizes from maxi- 
mally L = 150 at p = 1 to up to maximally L = 1000 at 
p = 0.52. At each pair of L and p we average typically 
well over about 10 4 instances. 

Before we proceed to collapsing the data, it is instruc- 
tive first to determine the exponent yp directly according 
to Eq. 10. Clearly, at p — p* , diluted lattices are almost 
all entirely reducible with the rules given in Sec. Ill Al and 
rarely is any subsequent optimization necessary. Hence, 
L is limited only by the cost of reduction itself, memory 
space, and statistics. The data for the defect energy at 
p* = 1/2 is plotted in Fig. 0] The asymptotic fit yields 
about yp — —0.98(2). This value seems to suggest that 
yp = —1, which would imply that the spin glass on a 
percolating cluster in two dimensions is essentially one- 
dimensional (where y = —1). Concluding that yp = — 1 
is exact may be misguided: Ref. 0] obtained yp ~ —0.99 
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FIG. 5: Plot on a logarithmic scale of the stiffness a as a 
function of system size L. The data is grouped into sets (con- 
nected by lines) parameterized by the bond density p. Most 
sets show a distinct scaling regime as indicated by Eq. for 
sufficiently large L. In particular, the data for p — 0.51 ap- 
pears to never reach scaling and will be disregarded entirely. 



FIG. 6: Scaling plot of the data from Fig. [K] for a, fitted 
to Eq. fllUP as a function of the scaling variable x — L/£(p) 
where = (p — p*)~ v ■ Data below the scaling regime in 
each set (i. e. below a certain value of L for each p) from 
Figs.|H|was cut. The continuous line represent the power-law 
[L/£(p)]"°- 281 , using the fitted value for y. 



on the basis of a scaling argument involving the numeri- 
cal solution of an integral. 

To obtain an optimal scaling collapse of the data in 
accordance with the discussion in Sec. IIIII we focus on 
the data in the asymptotic scaling regime for each set. 
To this end, we chose for each data set a lower cut in L 
by inspection of the data in Fig. [SJ All data points below 
the cut for each p are discarded, all data above are kept. 
Then the remaining data for all L and p are fitted to the 
four-parameter scaling form in Eq. (|l(Jfl . The resulting 
collapse is displayed in Fig. 

Initiating all of the parameters at near-optimal values, 
a least-square fit incorporating all the data shown in the 
collapse converges. The fitted values are y = —0.281(3), 
y P = -0.77(5), p* = 0.514(5), and /(oo) w 3.3. Errors in 
this fit are estimates based the sensitivity on varying each 
parameter and have to be judged cautiously. It should 
be noted how essential the inclusion of the parameter yp 
was for the collapse, even obtaining a fitted value not 
too far from its actual value determined in Fig. 0] The 
discrepancy between the fitted value yp = —0.77 and the 
accurate value yp = —0.99 is due to the unknown scaling 
for £(p) away from p c and the scaling-corrections for the 
approach Eq. J5J), as discussed above. 



of an advance than bond-diluted lattices did for higher- 
dimensional lattices. For one, in two dimensions there 
is no finite-temperature glass transition (y < 0) and 
conventional studies at full connectivity are successful 
at reaching large lattice sizes already, avoiding the ad- 
ditional uncertainties of a multi-parameter fit. In this 
scenario, Ref. |3lJ argued that such a collapse of data 
may provide diminishing returns for the computational 
effort. Similarly, it was observed there that a continu- 
ous bond distribution, considering the smaller size of ele- 
mentary excitations under bond-reversal, leads to larger 
scaling corrections in L. Those scaling corrections are 
enhanced further by open boundaries, which have been 
observed previously to result in only weakly decreasing 
corrections [71 l49l l50j . 

Our results indicate the validity of the recently pro- 
posed reduction scheme to determine the stiffness expo- 
nent y. Note that the fact that reduction works well in 
two-dimensions, where T c = holds, does not imply def- 
initely that it should work for d > 2. Nevertheless, since 
the overall behavior in d = 2 and higher dimensions is 
similar, it is highly probable that the reduction scheme is 
applicable also for higher dimensions |29l l30l | , where no 
exact ground-state algorithms are available. 



V. CONCLUSION 



We have used a scaling Ansatz proposed in Refs. |32j, 
|45| in conjunction with the spin reduction scheme of 
Refs. |23, Ho, H3 and exact ground-state calculations to 
study the defect energy at T = for bond-diluted lat- 
tices in two dimensions. The results for the stiffness ex- 
ponent y scale over 2 decades and are consistent with 
previous studies, validating the basic Ansatz. Yet, the 
obtained value y = —0.281(3) is only of comparable accu- 
racy to those studies, and the data collapse provides less 
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